Frustrated colloidal ordering and fully packed loops in arrays of optical traps 
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We propose that a system of colloidal particles interacting with a honeycomb array of optical 
traps that each contain three wells can be used to realize a fully packed loop model. One of the 
phases in this system can be mapped to Baxter's three-coloring problem, offering an easily accessible 
physical realization of this problem. As a function of temperature and interaction strength, we find 
a series of phases, including long range ordered loop or stripe states, stripes with sliding symmetries, 
random packed loop states, and disordered states in which the loops break apart. Our geometry 
could be constructed using ion trap arrays, BEG vortices in optical traps, or magnetic vortices in 
nanostructured superconductors. 

PACS numbers: 82.70.Dd,75.10.Hk 



Introduction. There has recently been tremendous 
growth in the area of creating idealized systems in which 
certain types of statistical mechanics models with and 
without geometric frustration can be physically realized, 
such as in nanomagnets and soft matter systems 

(3-[l4|. The key advantage of these systems is that they 
allow direct experimental access to the microscopic de- 
grees of freedom. One of the most active of these areas 
has been artificial spin ices created using nanomagnetic 
arrays with square fl^ or hexagonal ordering [H, Q , where 
ordered or frustrated states can occur that mimic real 
spin ice systems [15]. Here, various types of excitations 
such as monopoles can arise, and the dynamics can be 
studied under an external field [3] . There are many other 
statistical mechanics models that exhibit frustration ef- 
fects, including loop models such as the famous Baxter's 



three-coloring model [16|, where only very limited work 



has been performed on proposed physical realizations, 
all of which involve atomic degrees of freedom 



The nanomagnetic systems have certain constraints that 
make it very difficult to realize many other types of sta- 
tistical mechanics models of interest. 

In this letter we propose that a system of colloidal par- 
ticles interacting with optical trap arrays can be used to 
realize fully packed loop models, and show that one of 
the resulting phases can be mapped to the three-coloring 
model. Loop models have been applied to a wide variety 
of physical systems, ranging from polymer physics [21u22.] 
and turbulence [23] to optics [2J and magnetism [25|- 
[27! , and a physical realization of an idealized loop model 
would be a major step in this field. Golloidal systems 
interacting with periodic optical arrays have been exper- 
imentally realized 0, [HI, [si-fioi [28j and shown to exhibit 
novel types of orderings depending on the nature of the 
substrate 0, El, 0, [lH . Beyond these static states, it is also 
possible to study a variety of dynamical processes such 
as the motion of kinks and antikinks [lol, [HI- Highly 
tailored optical trap arrays can be created where the col- 
loidal particles can sit in multiple positions in a single 




FIG. 1: (a) Schematic diagram of the basic unit cell with two 
triple well traps each containing one colloidal particle, (b) 
and (c) are snapshots of a small portion of the system. The 
green triangles represent the traps and the red dots denote 
the particles, (b) shows a random distribution of particles 
at high temperatures. Uo pairs are circled and Ui pairs are 
boxed, (c) shows an example of a particle configuration that 
can be mapped to random fully packed loops in the hexagonal 
lattice, as illustrated in panel (d). The yellow contour in (d) 
corresponds to a flippable type-II loop. 



trapping site [H, HH, Hoj, so that arrays where colloidal 
particles can occupy one of three possible positions in a 
trap are well within current experimental capabilities. 

Model. We consider a 2D array of triangular-shaped 
traps that each contain three potential minima, as illus- 
trated in Fig. [IJa). These traps are similar to those ere- 



ated experimentally in Ref. Q . The traps form a honey- 
comb lattice with triangles of opposite orientations occu- 
pying the two different sublattices, as shown in Fig. 
and (c). When each trap contains one colloidal particle, 
the system provides a natural realization of the 3-state 
Potts model on the honeycomb lattice. We introduce a 
Potts variable <Ji = 1, 2, or 3 to denote the potential 
well occupied by the particle in the i-th trap. In order 
to characterize the colloidal ordering, we also introduce 
three unit vectors for each Potts state: ei = (0,1), and 
e2,3 = (±\/3/2, — 1/2). The particle in the i-th trap is 
located at = + d^^ , where is the center of the 
trap and the displacement d^r = i^^e^r, with the plus 
(minus) sign for up (down) triangles, and S denotes the 
linear size of the trap. 

The colloidal particles interact with each other via a 
repulsive screened Coulomb or Yukawa potential given by 
V{rij) = Voq'^ ex.-p{—K.rij)/rij. Here rij is the distance 
between two particles, Vq = Z*^/(47reeo), is the unit 
of charge, e is the solvent dielectric constant, q is the 
dimensionless colloid charge, and 1//^ is the screening 
length. The Hamiltonian of the model system reads 

n = ^^v{\R,i+d,^-d,^\) (1) 

where the summation runs over all pairs of triangular 
traps j. Since the particles are always confined to one of 
the three potential wells in each trap, we can identify the 
first few neighboring interaction terms of the colloidal po- 
tential V{rij) as summarized in Fig.O The dominant Uo 
interaction is between two particles at the closest corners 
of two neighboring triangles shown in Fig. Efa). Since 
these two potential wells are specified by the same Potts 
state in the respective traps, the Uq term essentially in- 
troduces an anisotropic antiferromagnetic interaction be- 
tween the Potts variables: 

= Uo^^Sa,j,a,^a,j,aj, (2) 
(ij) 

where (ij) denotes two nearest-neighbor traps, and aij = 
1,2,3 specifies the relevant Potts state of the adjacent 
wells of the (ij) pair. Such /7o-pairs of particles [circled 
in Fig. mb)] are energetically unfavorable and will be 
suppressed at temperatures T <C t/o- It is worth not- 
ing that the interactions in Eq. (j2j) are frustrated and 
there exist extensively degenerate Potts states (colloidal 
configurations without /7o-pairs) that minimize Hq. 

By attaching an arrow to each particle pointing from 
the center of the triangular trap to the corner occupied by 
the particle, the colloidal configuration can be mapped to 
a collection of directed strings. Since the triangles form a 
honeycomb lattice, a similar mapping can be established 
by extending the arrow onto the corresponding bond [see 
Fig. [TJc) and (d)]. As each trap contains exactly one 
particle, there is always an outgoing arrow for each ver- 
tex of the honeycomb lattice; however, the number of 
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FIG. 2: Various interaction terms arising from the screened- 
Coulomb or Yukawa potential V{rij) between a pair of colloids 
in the optical traps. 



incoming arrows for individual vertices can be 0, 1, or 
2. The number of vertices with no incoming arrow must 
equal the number of vertices with 2 incoming arrows since 
these are the sources and sinks (or end points) of the di- 
rected strings. The second and most relevant Ui term of 
the interaction, shown in Fig. EJb), prevents the fusion 
of two strings by penalizing vertices with two incoming 
arrows. Examples of Ui pairs are highlighted by square 
boxes in Fig. [IJb). The number of end point vertices is 
suppressed at temperatures T <C t/i, where for systems 
with periodic boundary conditions (BC) it becomes ener- 
getically more favorable for strings to form closed loops 
as shown in Figs. [TJc) and (d) [30|, [3l|. For finite lat- 
tices with open BC, the end points of the strings reside 
at the boundaries of the system. The further-neighbor 
interactions U2 and Us shown in Fig.[2](c,d) induce long- 
range ordering of particles at very low temperatures. In 
particular, the U2b term favors alignment of particles in 
two different alternating Potts states along one of the C3 
symmetry directions, effectively introducing a bending 
stiffness to the strings. 

It is worth noting that each fully packed loop (FPL) 
configuration on the honeycomb can be further mapped 
to a three-colored configuration on the same lattice. In 
Baxter's three coloring model [l6|, each bond of the hon- 
eycomb lattice is assigned a color R, G, or B, so that 
three different colors meet at each vertex, and all such 
configurations are given equal statistical weight. The R 
and B colored bonds thus form a FPL configuration as 
illustrated in Fig. md), and the two different sequences 
RBRB • • • and BRBR • • • correspond to the forward and 
backward propagating loops, respectively. It is impor- 
tant to note that all three-colored configurations are en- 
ergetically degenerate if we retain interactions up to the 
Ui terms only. Long-range orderings are induced by the 
further-distance interactions in V{rij). 

To understand the various phases of the model system, 
we perform Monte Carlo simulations on the effective 3- 
state Potts model described by Eq. ([T]). At high temper- 
atures, the standard single-site Metropolis updates are 
sufficient to equilibrate the system; however, such local 
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updates experience a dynamical freezing at temperatures 
T <C t/i due to the huge energy cost of updating a single- 
site Potts state. Since the effective degrees of freedom in 
this temperature regime are the fully packed loops dis- 
cussed above, we also implement two types of non-local 
updates in our Monte Carlo simulations similar to the 
loop alg orithm introduced for the three-coloring prob- 
lem [321 . In the first type of loop update, we randomly 
select a loop of head-to-tail arrows, or a RBRB • • • loop, 
and flip all the arrows; this move is accepted according to 
the standard Metropolis conditions with further-distance 
interactions 1/2^ Us^ • • • taken into account. The type-II 
loops consist of alternating bonds with and without ar- 
rows; they correspond to the GB or GR-colored loops in 
the 3-color scheme. An example type-II loop is shown in 
Fig.lHd). 

At very low temperatures, even the loop updates suffer 
freezing problems. Unlike the loops in dimer or spin-ice 



models [33|, [3^ , which can be constructed step by step 
from numerous possible paths, the loops in the three- 
coloring problem are predetermined by the colors in a 
given configuration. The so-called worm algorithm (sHj . 
in which detailed balance is always satisfied when con- 
structing the loop, cannot be applied to our case. The 
freezing problem arises because the acceptance rate of 
flipping a long loop in the standard Metropolis criterion 
becomes exceedingly small at low T. To overcome this 
problem, we employ the parallel tempering algorithm [36| 
to simulate the low-temperature regime. By simulta- 
neously simulating 150 replicas covering a temperature 
window < T < O.lt/i, we are able to fully equilibrate 
a system with periodic boundary conditions containing 
N = 2 X 6 X 12 particles; the results from a system with 
linear trap size (5 = 0.9 X (a/2V3) and screening length 
tv~^ = 0.06a, where a is the lattice constant of the un- 
derlying honeycomb lattice, are summarized in Fig. [3l 

In Fig. [3l^a) we plot the fraction Ny of honeycomb lat- 
tice vertices of type v as a function of temperature in 
the regime T <^ Uq. Since the occurrence of /7o-pairs is 
strongly suppressed in this regime, there exist only three 
vertex types Ndi^ Nd2^ and A^sc, defined according to the 
'coloring' of the three bonds attached to the vertex, as 
illustrated in the insets of Fig. [3l^a). The three bonds 
meeting at the lowest-energy N^^c vertices always have 
different colors. In the language of loops, these 3-color 
vertices have exactly one incoming and one outgoing ar- 
row. There are two types, A^^^i and A^d2, of higher-energy 
defect vertices that violate the three-color constraints; 
they correspond to the sources and sinks of the open 
strings, and always satisfy Ndi = Nd2' In Fig.[3l^a), as T 
decreases the fraction of defect vertices A'^^i + Nd2 grad- 
ually decreases before vanishing for T < Tsc ^ 0.1/7i, 
while the fraction of 3-color vertices Nsc saturates to 1 at 
low T. The system can thus be mapped to a three-colored 
or fully-packed loops configuration below the character- 
istic temperature T^c 
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FIG. 3: (a) Ny, the fraction of vertices of type v, as a func- 
tion of temperature T/Ui. Upper red line: Asc; lower blue 
line: A^^i + AcZ2; dashed line: all other vertex types. In- 
sets: schematics of the three low-temperature vertex types 
Aczi, Ad2, and Nsc. (b) Order parameters M and S along 
with Nsc as a function of temperature T/Ui. The parameter 
M characterizes a uniform long-range ordering of particles in 
which all loops are directed in the same direction and parallel 
to each other. The stripe order parameter S describes a par- 
tially ordered phase in which loops are parallel to each other 
but the direction of individual loops is disordered, (c) Phase 
diagram of temperature T in units of Voq^ /a vs showing 
the regions in which the ordered, smectic, three-coloring, and 
disordered states are observed. 



As discussed previously, the further-neighbor interac- 
tions U2 and Us induce long-range orderings of loops at 
lower temperatures. In particular, the loops acquire a 
bending stiffness due to the U2b interaction. As a result, 
the loops start to align themselves with one of the three 
principle lattice symmetry directions upon lowering the 
temperature. Since the dominant Potts interaction Uq 
is antiferromagnetic, we consider a Neel type order pa- 
rameter: M = {2/V3N)^.{-iye^^, where (-1)^ = +1 
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for up triangles and —1 for down triangles. The order 
parameter M = |M| indeed rises to its maximum at 
T <Tn ^ OmSUi as shown in Fig. [Sfb), indicating a 
ground state with long-range antiferro-Potts order. One 
of the perfectly ordered states is illustrated in Fig. HJ^a); 
there are a total of 6 degenerate ground states related to 
the breaking of Z2 sublattice (the arrows in the loops) 
and Cs rotational (the orientation of the loops) symme- 
tries. 

Interestingly, for decreasing temperature the order pa- 
rameter M shows an upturn at Ts ~ 0.034/7i, above 
the onset of long-range Potts order. Examination of 
the snapshots from Monte Carlo simulations shows a 
partially ordered phase with additional sliding symme- 
tries [37|. In this phase, the loops are either paral- 
lel or antiparallel to each other, hence breaking the 
Cs lattice rotational symmetry. The directions of ar- 
rows in individual loops remain disordered as shown in 
Fig. m^b). This partially ordered phase is characterized 
by a Z3 order parameter indicating the overall orienta- 
tion of loops and a set of Ising variables {tq, ti, • • • ^tl} 
specifying the direction of each loop. To characterize 
this stripe-like order, we first compute the antiferro- 
Potts order on a ID chain along one of the Cs axes: 
^a(c) = (l/L)vc^-^^^^(— l)^ecr^, where n is a site index 
along the chain c; a = 1, 2, 3 specifying the orientation of 
the chains; and Vq, = e/3 — e^, where {af3j) is a cyclic per- 
mutation of (123). The vector is used to project the 
vector sum to the relevant Potts states along the chain. 
Averaging over chains of the same orientation a gives a 
quasi- ID order parameter: = (1/^) l^al^)!^ 
finally the stripe order parameter is defined as their max- 
imum S = maxQ, Ma. As shown in Fig. [3l^b), the system 
enters the partially ordered stripe phase at T < as 
the order parameter S saturates to its maximum. 

We summarize the sequence of thermodynamic trans- 
formations, illustrated in the phase diagram in Fig. 3(c), 
as follows. As the temperature is lowered, the colloidal 
system first undergoes a crossover into the 3-color or ran- 
dom FPL phase at Tsc ^ 0{Ui). A phase transition 
into the partially ordered phase occurs at Ts ~ 0{U2h) 
when the stripe-ordering arises from the positive bend- 
ing energy produced by the U2h interaction. Finally, the 
system undergoes another phase transition into the long- 
range antiferro-Potts ordered ground state at Tn- We 
note that for larger system sizes, our Monte Carlo sim- 
ulations combining local Metropolis, loop updates, and 
parallel tempering are able to reach the equilibrium 3- 
color phase at T < Tsc However, since full equilibration 
to the partially ordered striped phase as well as the fully 
ordered ground state requires fiipping system-size loops, 
which costs too much energy for larger lattices, our al- 
gorithm can only produce a multi-domain stripe phase. 
It is worth noting that in the thermodynamic limit, the 
system cannot reach the true long-range order and stays 
in this smectic-like phase due to the huge energy barrier 




(a) (b) 

FIG. 4: (a) A long-range ordered loop state characterized by 
an antiferromagnetic Potts order parameter M. The parallel 
loops in this ordered state are directed in the same direction, 
(b) A partially ordered loop state exhibiting a sliding symme- 
try. The loops in this phase are parallel to each other but the 
direction of individual loops remains disordered. This state 
is characterized by the stripe order parameter S. 



separating different stripe states. 

In summary we have proposed that colloidal particles 
interacting with a honeycomb array of optical traps that 
each contain three wells can be used to realize a fully 
packed loop model. We show that this system exhibits 
an ordered ground state, a smectic-like stripe phase with 
a sliding symmetry, a random fully packed loop state, 
and a disordered state with broken loops. The random 
fully packed loop state can be mapped to Baxter's three- 
coloring problem, indicating that our system could be 
used to create a physical realization of this problem. We 
map out where these phases occur as a function of tem- 
perature and interaction strength. Fully packed loops 
on different lattices can be similarly realized with optical 
arrays in which the number of potential wells in a trap 
site is the same as its coordination number. Our results 
should be generalizible to other systems of repulsively in- 
teracting particles in a similar array of three-well traps, 
such as for vortices in EEC's interacting with optical ar- 
rays, vortices in nanostructured type-II superconductors, 
and ions in tailored trap arrays. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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